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CROSSTEX  -  Wave  breaking,  boundary  layer  processes,  the  resulting 
sediment  transport  and  beach  profile  evolution 

Tian-Jian  Hsu' 

Civil  and  Coastal  Engineering,  University  of  Florida 
Applied  Ocean  Physics  and  Engineering,  Woods  Hole  Oceanographic  Institution^ 

Award  Number:  N00014-07-1-0154,  NOOO 14-05- 1-0082 

LONG-TERM  GOALS 

To  develop  and  test  with  laboratory  and  field  data  a  robust  modeling  framework  that 
predicts  hydrodynamics  and  the  fate  of  terrestrial  and  marine  sediment  in  the 
heterogeneous  environment. 

OBJECTIVES 

■  Identify  the  role  of  steep  waves  and  breaking  wave  turbulence  on  sediment 
suspension  events  and  its  implication  on  on/offshore  sediment  transport  rate. 

■  Develop  phase/depth-resolving  numerical  models  for  bottom  sediment  transport  and 
surf  zone  hydrodynamics.  Validate  these  models  using  data  measured  in  CROSSTEX 
and  other  experiments. 

■  Develop  simplified  phase-resolving  formulations  for  concentrated  sediment  transport, 
suspended  load  transport  and  its  near-bed  boundary  conditions  under  breaking  waves. 

APPROACH 

A  two-phase  flow  model  (Hsu  et  al  2004;  Amoudry  et  al.  2008)  is  extended  to  model 
sand  transport  driven  by  measured  random  wave-current  forcing  during  CROSSTEX 
(Scott  et  al.  2007).  Measured  near-bed  time  series  of  flow  velocities  are  analyzed  to 
extract  (approximately)  the  time  series  of  wave  velocities  and  turbulence  quantities  (TKE 
and  dissipation  rate).  These  time  series  are  then  used  to  drive  the  two-phase  flow  model 
and  to  calculate  and  resulting  sand  transport  from  immobile  sand  bed  to  the  dilute  region 
of  about  5  cm  above  the  immobile  bed.  Physical  experiments  also  provide  measured 
suspended  sand  concentration  in  the  intermediate  to  dilute  region  (volume  concentration 
~<15%),  which  can  be  used  to  identify  the  suspension  event  and  to  validate  the  numerical 
model.  On  the  other  hand,  two-phase  model  results  for  concentrated  sediment  dynamics 
can  bridge  the  missing  information  at  the  regime  of  concentration  from  about  15%  to 
random-close  packing  (~63%),  which  is  difficult  to  measure  in  the  physical  experiment. 

Because  it  is  non-trivial  to  separate  turbulence  from  the  wave  motion  in  the 
measured  velocity  time  series,  also  because  it  is  perhaps  unclear  at  this  point  what  is  the 
major  cause  (e.g.,  steep  wave  or  breaking  wave  turbulence)  of  the  strong  sediment 
suspension  event,  it  is  critical  to  utilize  a  more  powerful  data  analysis  tool  than 
commonly  used  Fourier  analysis.  Wavelet  analysis  is  able  to  identify  not  only  the  steep 
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oscillation  events  but  also  their  (temporal)  scale  (Scott  2005;  Scott  et  al.  2005).  Hence, 
more  information  can  be  extracted  from  the  measured  time-series  when  using  wavelet 
analysis  in  conjunction  with  Fourier  analysis.  Wavelet  analysis  on  the  CROSSTEX  data 
is  investigated  by  Research  Scientist  Dr.  Nicolas  Scott  of  University  of  Florida  funded  by 
this  project. 

WORK  COMPLETED 

CROSSTEX  is  a  collaborative  research  project  aiming  at  quantifying  various 
mechanisms  of  wave-  and  current-induced  sediment  transport  in  the  surf  zone.  The 
physical  experiment  of  CROSSTEX  was  completed  in  Sep  2005  at  O.  H.  Hinsdale  Wave 
Laboratory  of  Oregon  State  University  (Scott  2006;  Scott  et  al.  2007).  The  experimental 
design  of  CROSSTEX  is  to  mimic  the  field  condition  of  onshore/offshore  sandbar 
migration  events  observed  during  Duck  94  (Gallagher  et  al.  1998)  so  that  possible 
mechanisms  controlling  the  onshore/offshore  sediment  flux  can  be  carefully  studied. 

Even  in  a  well-controlled  laboratory  facility,  detailed  sediment  transport  in  the 
concentrated  regime  near  the  bed  cannot  be  measured  with  high  temporal  and  spatial 
resolution.  In  this  study,  we  utilize  a  two-phase  model,  which  resolves  both  the  bedload 
and  suspended  load,  to  link  the  missing  information  in  the  concentrated  regime  of 
transport.  The  two-phase  model  for  sand  transport  has  been  validated  in  the  past  two 
years  with  published  laboratory  U-tube  results  (Amoudry  et  al.  2008;  Hsu  and  Yu  2008). 
This  year,  we  also  extend  the  model  with  nonlinear  boundary  layer  terms  (i.e., 
approximate  djdx  »  c"*  djdt  in  the  2DV  two-phase  equations  with  c  the  wave  celerity)  to 
mimic  sediment  transport  under  progressive  waves.  The  model  is  validated  with 
SISTEX99  data  (Dohmen- Janssen  and  Hanes  2004)  in  order  to  investigate  the  effect  of 
nonlinear  boundary  streaming  on  wave-induced  transport  (Trowbridge  and  Young  1989; 
Henderson  et  al.  2004)  in  comparison  with  wave  skewness/asymmetry  (e.g.,  Hsu  and 
Hanes  2004).  It  is  found  that  for  the  grain  diameter  (0.24  mm)  and  wave  condition  used 
during  SISTEX99,  nonlinear  wave  shape  (e.g.,  skewness)  is  the  dominant  mechanism  for 
onshore  transport;  nevertheless,  nonlinear  boundary  layer  streaming  effect  causes 
additional  36%  of  the  onshore  transport.  Numerical  experiments  further  suggest  that 
nonlinear  streaming  effect  becomes  increasingly  important  for  finer  grain  and  nonlinear 
boundary  layer  process  can  increase  onshore  transport  in  the  shoaling  zone  and  bar  crest 
regime  by  about  80~90%,  which  contributes  a  non-negligible  portion  of  onshore  transport 
to  that  purely  due  to  nonlinear  wave  shape.  Research  results  on  this  topic  are  summarized 
in  a  manuscript  that  shall  be  submitted  for  publication  in  the  near  future. 

Starting  in  FY08,  CROSSTEX  data  is  analyzed  using  wavelet  transform  and 
statistical  analysis.  Morlet  wavelet  transform  has  been  previously  shown  to  be  effective 
in  identifying  steep  wave  in  open  ocean  (Scott  et  al.  2005).  Similar  method  is  applied 
here  to  the  measured  velocity  time-series  (50  Hz  sampling  rate).  Measured  cross-shore 
(x)  velocity  component  has  both  wave  and  turbulence  information  in  the  frequency  range 
1-2  Hz,  while  the  alongshore  (y)  velocity  component  has  clean  turbulence  information. 
This  suggests  that  the  wavelet  analysis  must  be  applied  concurrently  to  x  and  y 
components  to  obtain  a  better  understanding  of  steep  wave  and  breaking  wave 
turbulence.  When  plotting  the  wavelet  transform  results  for  velocities  along  with 
measured  sediment  concentration  time  series  at  given  location  above  the  bed,  we  can 


identify  suspension  events  correlated  with  different  scenarios  of  steep  wave  and/or 
breaking  wave  turbulence  events  (see  Figure  1).  Our  goal  is  to  come  up  with  useful 
statistics  on  the  occurrence  frequencies  of  sand  suspension  events  due  these  scenarios 
(see  Table  1).  This  information  can  be  used  to  interpret  model  results  calculated  by  the 
1 D V  two-phase  model  that  only  captures  locally  generated  events  as  opposed  to  nonlocal 
advective  processes.  Research  results  on  this  topic  were  presented  in  2008  AGU  Ocean 
Science  Meeting  and  31*'  International  Conference  on  Coastal  Engineering.  In  addition,  a 
manuscript  is  also  submitted  for  publication. 

RESULTS 

Figure  la  shows  a  U  component  velocity  time  series  from  the  erosive  case,  a  contour  plot 
of  its  wavelet  transform,  and  the  concurrently  measured  sediment  concentration  time 
series.  The  asterisks  identify  high  acceleration  events  at  small  and  large  frequencies  that 
have  accelerations  above  0.3  m/s^  and  they  are  visibly  correlated  with  the  local  maxima 
of  velocity  variation  in  the  velocity  signal.  The  algorithm  is  able  to  perform  a  multi-scale 
filtration  of  the  data  to  reveal  high  acceleration  events  associated  with  cohesive  wave 
groups  at  both  large  and  small  scales.  Figure  lb  shows  a  V  component  velocity  time 
series  from  the  erosive  case,  a  contour  plot  of  the  wavelet  transform,  and  the  sediment 
concentration  time  series  taken  from  the  same  instant.  The  asterisks  identify  high 
turbulent  velocity  events  that  have  a  velocity  above  0.03 m/s .  Combining  results  in  both 
Figure  la  and  lb,  we  can  identify  a  typical  scenario  that  a  steep  wave  causing  significant 
breaking  wave  turbulence  approaching  the  bed  and  induces  strong  sand  suspension. 

Table  1  shows  the  steep  wave,  turbulence,  and  concentration  statistics  for  the 
erosive  and  accretive  conditions  at  1  cmab  obtained  using  the  wave-breaking  advective 
timescale,  which  is  taken  to  be  the  ratio  of  average  breaking  wave  height  and  RMS  wave 
velocity  associated  with  each  case.  The  top  part  of  Table  1  show  the  statistics  of  the 
scenarios  (shown  in  the  1*'  column)  associated  with  the  passage  of  steep  waves.  On  the 
other  hand,  the  bottom  part  of  Table  1  comprises  the  statistics  of  the  scenarios  associated 
with  the  occurrence  of  high  sediment  concentration  events.  During  the  erosive  condition, 
only  3 1  %  of  the  steep  wave  events  have  breaking  wave  turbulent  events  that  approach  1 
cmab  (see  the  3'^  row  of  Table  1).  For  accretive  case,  this  statistic  further  decreases  to 
18%.  In  addition,  for  the  erosive  condition,  only  14  %  of  steep  waves  are  associated  with 
turbulence  and  high  concentration  events  while  this  statistic  is  5  %  for  accretive 
condition.  This  last  statistical  category  shows  that  based  on  local  point  measurements 
within  the  wave  boundary  layer  for  both  cases,  the  percentage  of  steep  waves  associated 
with  turbulence  and  high  concentration  events  is  a  small  fraction  of  the  total  amount  of 
steep  waves.  The  classic  scenario  that  steep  breaking  waves  cause  high  turbulence 
approaching  the  bed  and  directly  induced  sand  suspension  is  also  a  fraction  of  the  total 
amount  of  sediment  suspension  events.  For  the  erosive  condition  30%  of  high 
concentration  events  are  associated  with  steep  wave  events  and  turbulent  events  (see  the 
7"'  row).  This  same  statistic  for  the  accretive  condition  drops  to  15%.  All  three  categories 
discussed  here  suggest  that  steep  waves  are  breaking  intermittently  causing  intermittent 
bottom  wave  boundary  layer  sediment  suspension.  The  intermittency  of  wave  breaking 
and  sediment  suspension  is  qualitatively  consistent  with  results  reported  in  the  literature 
(e.g.  Cox  and  Kobayashi  2000).  It  is  also  noted  that  while  the  accretive  case  has  a  larger 


number  (33)  of  steep  waves  exceeding  acceleration  threshold  than  the  erosive  case  (27.7), 
the  erosive  case  has  a  larger  percentage  of  steep  waves  associated  with  high 
concentration  and  turbulence  events.  In  other  words,  erosive  case  appears  to  have  more 
intense  or  localized  (the  sensors  are  collocated  at  one  cross-shore  location  x=63.7  m) 
breaking  wave  turbulence  and  more  frequent  suspension  events  than  those  of  accretive 
case  despite  the  fact  that  the  RMS  wave  velocity  is  greater  for  the  accretive  case.  Further 
analyses  suggest  that  these  statistics  are  sensitive  to  the  choice  of  time  window  size  for 
accretive  condition  and  measurement  location  above  the  bed  (for  both  erosive  and 
accretive  condition).  It  is  concluded  in  this  study  that  when  large  time  window  size  is 
utilized  and  data  measured  at  higher  locations  are  studied,  more  nonlocal  processes,  such 
as  turbulence  and  sediment  advected  from  the  upstream,  can  change  the  results. 
Therefore,  statistical  results  shown  in  Table  1  are  thought  to  better  represent  locally- 
generated  processes,  which  is  consistent  with  the  IDV  two-phase  modeling  results.  It  is 
also  note  here  that  according  to  Table  1,  only  about  1/3  of  the  suspension  events  are 
locally-generated  and  significant  portion  of  the  suspension  events  are  possibly  due  to 
nonlocal  advective  processes. 

Numerical  modeling  results  for  total  load  transport  display  agreement  with 
observed  laboratory  measurements  of  bathymetry  evolution  for  the  erosive  condition.  The 
model  incorporating  breaking  wave  turbulence  is  able  to  reproduce  suspension  events 
caused  by  local  steep  wave  and  breaking  wave  turbulence  that  give  sand  concentration 
profiles  in  agreement  with  measured  data.  Model  results  suggest  mean  undertow  current 
causes  offshore  transport  which  work  against  the  onshore  directed  wave-induced 
transport  to  finally  give  a  net  offshore  sediment  flux.  Model  results  also  suggest  under  a 
steep  wave  breaking  event,  the  breaking  wave  turbulence  can  enhance  net  transport  by 
50-150%.  However,  the  direction  of  the  enhanced  transport  is  determined  by  the  relative 
timing  between  the  breaking  wave  turbulence  event  and  wave  phase.  Correlation  and 
statistical  analysis  of  measured  data  suggest  that  about  2/3  of  the  high  breaking  wave 
turbulence  events  are  negatively  correlated  with  the  wave  velocity  (see  the  last  row  of 
Table  1).  Hence,  our  data  analysis  support  breaking  wave  turbulence  tends  to  enhance 
offshore  sand  transport.  Model  results  of  net  transport  rate  disagree  with  the  measured 
(small)  beach  profile  change  for  accretive  condition.  The  onshore  sediment  transport  in 
present  IDV  numerical  model  is  mainly  caused  by  nonlinear  wave  shape,  such  as  wave 
velocity  skewness  and  asymmetry.  It  is  possible  that  other  onshore  sediment  transport 
mechanism,  such  as  non-local  advection  and  nonlinear  boundary  layer  streaming,  which 
are  not  captured  by  the  present  model,  are  likely  mechanisms  for  the  onshore  sediment 
transport  observed  in  this  wave  flume  experiment  during  accretive  condition. 

IMPACT/APPLICATIONS 

The  present  research  efforts  focus  on  new  data  analysis  techniques  for  surf  zone  sediment 
transport  developing/validating  detailed  numerical  models  for  sediment  transport  and 
wave  hydrodynamics.  We  have  focused  on  refining  the  numerical  schemes  in  the  IDV 
two-phase  model  and  2DV  RANS  wave  model  so  that  they  will  be  more  robust  and  user- 
friendly  in  the  future.  Recent  effort  using  wavelet  analysis  to  study  sediment  transport 
under  breaking  wave  is  also  new  to  the  coastal  sediment  transport  community. 
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Figure  1:  Wavelet  transform  for  U  (represents  wave)  and  V  (represent  turbulent  fluctuation) 
velocity  components  at  1  cmab  during  erosive  condition,  (a)  Wavelet  transform  for  U  component 
with  asterisks  delineate  high  acceleration  events.  Acceleration  threshold  of  0.3  m/s^  is 
designated  by  dotted  line.  Contours  drawn  at  every  0.12  interval,  (b)  Wavelet  transform  for  V 
component  with  Asterisks  delineate  high  velocity  events.  Velocity  threshold  of  0.03  m/s  is 
designated  by  dotted  line.  Contours  drawn  at  every  0.017  interval.  Thick-solid  curve  in  the 
wavelet  transform  plots  are  corresponding  measured  concentration  time  series  at  Icmab. 


Erosive  Case 

Runl 

Run2 

Run3 

Average 

Accretive  Case 

Total  #  of  Steep 

28.6 

27.2 

27.4 

27.7  ±1.4 

Wave  Events 

33.75 

32.2 

33  ±3.9 

%  of  steep 

8.6,  30  % 

7.4,  27% 

10.2,  37% 

31  %  ±  7.0% 

waves  with 
turbulent  events 

6.25,19% 

5.4,  17% 

18%  ±6.2% 

%  of  steep 

6.8,  24% 

7,  26% 

7,  26%. 

25  %±8.7% 

waves  &  high 
cone,  events 

4.25,13% 

3.4,  11  %. 

12%o±5.1%o 

%  of  steep 

3.8,  13% 

3.4,  13% 

4,  15% 

14%  ±4.6% 

waves  with 
turbulent  &  high 
cone,  events 

2.5,  7% 

1.2,  3.7% 

5%±3.6% 

Total  #  of  High 

16.2 

19 

18.2 

17.8±6.8 

Cone.  Events 

12.75 

8.6 

10±4.0 

%  of  high  cone. 

5.2,  32% 

4.8,  25%. 

6.0,  33%. 

30%±10%o 

events  with 
steep  & 

2.75,22% 

0.6,  7% 

15%  ±6.2%. 

turbulent  events 

%  of  high  cone. 

3.2,  20% 

4.8,  25% 

4.2,  23  % 

22  %  ±  6.9% 

events  near 
nothing 

4.75,  37% 

4.2,  49% 

43%  ±15% 

%  of  high  cone. 

4.8,  30  % 

4.6,  24  % 

4.6,  25% 

26%  ±9.1% 

events  near  steep 
waves  only 

3.5,  27% 

1.6,  19% 

23%.±12%o 

%  of  high  cone. 

3,  19% 

4.8,  26  % 

8.8,  19% 

21  %  ±13.6%o 

events  near 

turbulent  events 

1.75,  13%. 

2.2,  26% 

19%  ±9.1% 

only 

%  of  turbulent 

68  %o 

74%o 

75%. 

72%o±7.5% 

events  in 
negative  phase 

65% 

61%. 

63%o±6.1% 

Table  1:  Statistics  of  each  scenario  for  steep  wave,  high  turbulence,  and  high  sand  concentration 
events  for  both  erosive  (upper)  and  accretive  (lower)  conditions.  Scenario  descriptions  are  shown 
at  the  1st  column  followed  by  the  resulting  statistics  for  each  run  (2"^  to  4*^  columns)  and 
averages  over  all  runs  (5^  column).  Statistics  are  for  time  series  measured  1  cmab.  Wave¬ 
breaking  advective  timescale  are  used  for  detection  window  size  (1.8  sec  for  erosive  case  and  1.2 
sec  for  accretive  case). 
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